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Gravitational waves emitted by binary systems in the inspiral phase carry a complicated structure, 
consisting in a superposition of different harmonics of the orbital frequency, the amplitude of each of 
them taking the form of a post-Newtonian series. In addition to that, spinning binaries experience 
spin-orbit and spin-spin couplings which induce a precession of the orbital angular momentum and 
of the individual spins. With one exception, previous analyses of the measurement accuracy of gravi- 
tational wave experiments for comparable-mass binary systems have neglected either spin-precession 
effects or subdominant harmonics and amplitude modulations. Here we give the first explicit de- 
scription of how these effects combine to improve parameter estimation. We consider supermassive 
black hole binaries as expected to be observed with the planned space-based interferometer LISA, 
and study the measurement accuracy for several astrophysically interesting parameters obtainable 
taking into account the full 2PN waveform for spinning bodies, as well as spin-precession effects. 
We find that for binaries with a total mass in the range 10 5 A'/ Q < M < 10 7 M Q at a redshift of 
1, a factor ~ 1.5 is in general gained in accuracy, with the notable exception of the determination 
of the individual masses in equal-mass systems, for which a factor ~ 5 can be gained. We also 
find, as could be expected, that using the full waveform helps increasing the upper mass limit for 
detection, which can be as high as M — 10 s Mq at a redshift of 1, as well as the redshift limit 
where some information can be extracted from a system, which is roug hly z > 10 for M < 10 7 M Q , 
1.5-5 times higher than with the restricted waveform. We computed that the full waveform allows 
to use supermassive black hole binaries as standard sirens up to a redshift of z » 1.6, about 0.4 
larger than what previous studies allowed. We found that for lower unequal-mass binary systems, 
the measurement accuracy is not as drastically improved as for other systems. This suggests that 
for these systems, adding parameters such as eccentricity or alternative gravity parameters could 
be achieved without much loss in the accuracy. 



I. INTRODUCTION 

Gravitational waves (GW's), once their observation is 
made possible, will provide new means of observing the 
Universe. So far, the vast majority of observations have 
been made through electromagnetic radiation, and grav- 
itational waves will surely make visible different aspects 
of the Universe. For example, one could probe the differ- 
ent galaxy formation models by detecting supermassive 
black hole mergers in a large redshift range [jj. Or, as 
GW's provide a good way of measuring the luminosity 
distance to their source, one could combine gravitational 
and electromagnetic observations to build a robust mea- 
surement of the Hubble diagram, which would be of great 
interest for cosmology 0, Q ■ Another possibility would 
be to measure alternative gravity parameters @, H> @- 
This can potentially be a powerful way to constrain such 
theories, as each observed GW will give an independent 
measurement of their parameters. 

The new generation of ground-based detectors, such 
as advanced LIGO, and the space-based detector LISA 
will probably make the direct detection of gravitational 
waves possible. Some of the most important sources of 
gravitational waves are the compact binary systems, i.e. 
systems of two compact objects (white dwarfs, neutron 
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stars, or black holes). As such detections rely on matched 
filtering techniques, several groups have made efforts in 
building accurate templates based on the post-Newtonian 
(PN) approximation (see e.g. @, i, ©, M, El El ) ■ The 
limitations of such results have then been checked to 
estimate the precision with which one could measure 
the properties of a system emitting such waves, for ex- 
ample its distance from the Solar System, location in 
the sky, or the individual masses of the objects forming 

it 0, EI El El El El El, HI- 

A binary system of compact objects emits gravita- 
tional waves during three distinct phases, called inspiral, 
merger, and ringdown. During the inspiral phase, most of 
the gravitational radiation is emitted at twice the orbital 
frequency of the system, which slowly increases as it loses 
energy emitting gravitational waves. As its two members 
come closer, higher harmonics become more and more im- 
portant, and more power gets emitted. Finally, the two 
members enter the merger phase, where they have come 
so close that they cannot be treated anymore as two sep- 
arate objects, and begin to merge, emitting complicated 
gravitational radiation that has not yet been described 
but numerically. After that, the remnant begins to ra- 
diate away its energy during the ringdown phase, where 
it approaches exponentially the structure of a Kerr black 
hole. 

LISA has been designed to be particularly sensitive to 
binaries that contain a supermassive black hole (SMBH, 
of mass 10 5 - 10 8 Af Q ), which can be separated into two 
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(or three) categories, according to the companion mass: 
extreme mass-ratio inspirals (EMRIs) and supermassive 
black hole binaries (SMBHBs) . A third category is some- 
times added: intermediate mass ratio inspirals (IMRIs), 
which somehow lies between the two others. 

EMRIs are inspirals with a mass ratio between the two 
members of the order of 10~ 4 - 10~ 7 , and are most accu- 
rately described by black hole perturbation theory. Such 
events are likely to occur in the center of galaxies, the 
majority of which are believed to host a SMBH, when a 
compact object of stellar mass is "eaten up" by the cen- 
tral object. Several groups have attacked the problem of 
describing the form of the gravitational radiation emitted 
by such objects [2ll. [22I [23j . 

SMBHBs are binaries containing two SMBHs, forming 
during the merging of two galaxies, when they both host 
one in their center. Such events should be much more 
rare than EMRIs, but some galaxy formation models [24| 
and observation of nearby galaxies [25[ suggest that they 
might happen often enough to be observable. As these 
events are much louder than any other source, we could 
observe them at very high redshifts (up to z = 20 or 
even higher) , and thus tightly constrain galaxy formation 
models. 

The first attempt to estimate with which accuracy an 
interferometer could measure the properties of a compact 
object binary was made by Finn [13| |. who first intro- 
duced the Fisher matrix analysis, which is now widely 
used in this context. A few years later, Cutler [l4| ap- 
plied this formalism to LISA, focusing on the angular 
resolution that the space-based detector could get for 
black hole binaries, using the Newtonian quadrupole for- 
mula. Hughes repeated the study including the PN 
expansion for the frequency of the wave. Vecchio flol . 
then, considered the case of the "simple precession" [26| 
of the angular momenta for spinning BH's. Lang and 
Hughes [17| then used the full precession equations to 
further refine the parameter estimation. Recently, Arun 
et al. [18| and Porter and Cornish [19( included the full 
post-Newtonian waveform in the context of nonspinning 
black holes, and Trias and Sintes [2(| used it for spin- 
ning black holes neglecting spin-precession effects. The 
LISA Parameter Estimation Taskforce [27| used the full 
waveform with spin-precession effects, without publish- 
ing a detailed study of the expected statistical errors. It 
is worth noting that all of the effects that these works 
studied, as more and more precise waveforms were used, 
helped to improve subsequently the expected measure- 
ment accuracy of LISA. In this paper, we study the ques- 
tion of whether the inclusion of the full post-Newtonian 
waveform in the context of spinning black holes undergo- 
ing spin-orbit precession helps to break more degenera- 
cies, thus helping to further increase the accuracy in the 
measurement of the source parameters. 

In Sec. UH we derive the waveform used to perform 
our study, and we quote some basics of Fisher matrix 
analysis. In Sec. IIII1 we describe the simulations that we 
ran. In Sec. IIV1 we give the results of our simulations, 



and analyze them from an astrophysical point of view. 
We conclude in Sec. |Vl 



II. THEORY 
A. Evolution 

The state of a binary system of two Kerr black holes 
at a given time in the center of mass frame is fully de- 
scribed by 14 intrinsic parameters. These reduce to 12 if 
we assume that the binary lies on a circular orbit. One 
possible choice is to take as intrinsic parameters a unit 
vector pointing in the direction of the orbital angular 
momentum £, the orbital angular frequency ui (we will 
reserve the symbol / for arguments of Fourier transforms, 
and will express the orbital frequency always with the an- 
gular frequency to avoid confusions), the individual spins 
of each black hole, Si and S^, their masses mi and 777,2, 
and the orbital phase (p. To these intrinsic parameters, 
we have to add three more extrinsic parameters which 
locate the binary in space. Those can be chosen to be 
n, a unit vector pointing in the direction of the binary 
as seen from the Solar System, and c?l, the luminosity 
distance from the binary to the Sun. We will denote all 
unit vectors with a hat throughout this paper. 

Another extrinsic parameter, the redshift, also plays 
a role in the determination of the waveform, but it can- 
not be detected by GW observations. Indeed, the red- 
shift causes the observed angular frequency f of the 
wave to decrease with respect to the emitted one f e , 
as fo — Je/(1 + z). But (see following derivation) 
the exact same wave, within the post-Newtonian frame- 
work, is emitted by a second system with parameters 
mf'' = (1 + z)rrSp , = (1 + z)d^\ not experiencing 
any redshift. Therefore, the redshift and luminosity dis- 
tance cannot be measured separately with a gravitational 
wave observation, so that we have to assume a relation 
between the two parameters. This implies that obser- 
vations of a light signal emitted during a merger, the 
redshift of which is possible to determine, are of great 
astrophysical interest. During the whole derivation of 
the gravitational wave signal below, we assume that the 
source is at redshift z = 0. The actual observed wave can 
then be easily determined redshifting the masses and lu- 
minosity distance, as we did in our simulations. 

To compute the relation between redshift and luminos- 
ity distance, we assume a flat ACDM cosmology with- 
out radiation with the latest WMAP parameters [28| : 
fl A = 0.72, tt m = 0.28, H = 70.1 km/s/Mpc. 

The relation is then given by 

c r z dz' 
d L (z) = (l + z)— , (1) 

Ho Jo y/ajY+z 7 ) 3 + n A 

which can be determined numerically. 

The problem of the motion of the system during the in- 
spiral phase in full General Relativity has been too hard 
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to be solved so far. However, a great effort has been 
made to attack the problem in the framework of the post- 
Newtonian formalism. The current state-of-the-art evo- 
lution equations go up to the 2.5PN order beyond leading 
order for spinning objects [9(. As 2.5PN spin-orbit and 
spin-spin coupling terms are not yet known in the wave- 
form, we chose to stop at the 2PN level, up to which both 
the evolution equations and the waveform are known. We 
will use the following mass parameters for the derivation 
of the evolution equations and of the waveform: the total 
mass M = mi + m 2 , the reduced mass fj, = mim^/M, 
and the symmetric mass ratio v = fi/M. 

The 2PN orbit-averaged relation between the orbital 
angular frequency uj and the orbital separation in har- 
monic coordinates r is given by fiol ] 
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where the orbital separation parameter 7 and the spin- 
orbit and spin-spin couplings j3 and a are given by 
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The evolution equation of the angular frequency is 
given at 2PN order by [10( 
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where x is a dimensionless orbital frequency parameter 
defined as 
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We can integrate Eq. ([6]) to get 
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Integrating once more yields the orbital phase if = 
J udt, as a function of the orbital frequency parameter 
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The dragging of inertial frames induces a coupling be- 
tween the individual spins and the orbital angular mo- 
mentum. The orbit-averaged conservative part of the 
evolution equations (i.e. without radiation reaction, 
L + Si + S 2 = 0) are given for circular orbits at 2PN 
order by 
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where it is understood that i ^ j, i,j S {1,2}, and the 
orbital separation r and the norm of the orbital angular 
momentum L are related to the orbital frequency by their 
Newtonian relation: 
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Higher order relations would give corrections which ex- 
ceed the 2PN order. 

Using the above relations together with the first order 
of Eq. ([H]), we can change variables from time to orbital 
angular frequency, and use the relations to express the 
precession equations: 
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where 



1. Extrinsic effects 
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Waveform 
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The general form of a gravitational wave emitted by 
a two-body system, even nonspinning, is not known in 
the context of full general relativity. However, it has 
been computed in the post-Newtonian framework. The 
results for a nonspinning binary system are available at 
2.5PN order [IH, and the spin effects at 2PN order [13]. 

A convenient way to define the phase of the wave ob- 
served in a detector is in terms of the "principal-l- direc- 
tion" |26[ , which is defined as the direction of the vector 
L x fi. As the orbital angular momentum precesses, the 
principal+ direction changes, and this must be taken into 
account in the waveform. This effect amounts, at 2PN 
order, to 
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where u>q is an arbitrary constant corresponding to the 
time to, Sipo = — f t "(dS(p/dt)dt, and dL/du> is given in 
Eq. dUD. 

The 2PN accurate orbital phase is then given in terms 
of orbital angular frequency by: (f)(u>) — (f(u>) + S(p(u). 

The waveform is a series of harmonics of the orbital 
frequency: 
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The coefficients of the series take the form of post- 
Newtonian series: 
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The exact form of the coefficients for a nonspinning 
system can be found in [ill, Note, however, that 

both express their final result using another phase which 
differs from the orbital phase at 1.5PN order: \t = (f> — 
2 \og(Lu/Cu)x 3 / 2 , where u> is an arbitrary constant. We 
put together the results from these two papers to build 
a coherent 2PN accurate waveform for spinning bodies, 
see the appendix. 



The LISA constellation will consist in three spacecrafts 
launched in orbit around the Sun, at a mean distance of 
1 AU, on slightly eccentric orbit so that the spacecrafts 
stay at the same distance from each other all along the 
year. The barycenter of LISA will be located on the orbit 
of the Earth, 20° behind it, and the normal to the plane 
on which the spacecrafts lie will make a 60° angle with 
the normal to the ecliptic, see Fig. [1] 

To describe extrinsic effects that depend on the po- 
sition of LISA, we follow and define two different 
frames: a frame tied to the detector, (x, y, z), and a fixed, 
Solar System frame, tied to the distant stars (x, y, z) (we 
consider that the motion of the Sun with respect to the 
distant stars can be neglected during the lifetime of the 
LISA mission). 

The unit vectors along the arms of LISA Zj, i E {1,2,3} 
are defined in the detector frame: 



cos 7^ + smj t y, 
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The (x,y) plane of the Solar System frame is defined 
to be the ecliptic, so that the spherical angles of the 
barycenter of LISA are 



= 2irt/T, 



(23) 



where T = 1 yr, and we chose that <!> = at t = 0. 

The waveform is given relative to the Solar System 
frame. To take into account the fact that the detector 
is not static in this frame, we have to add a phase to 
each harmonic, which is equivalent to add the so-called 
Doppler phase to the orbital phase: 



(t>D{t) = sm0 N cos(^>(t) -fa), (24) 

c 



where R — 1 AU, and 6m and 4>n are the spherical angles 
of the position of the source in the Solar System frame. 

The orbital phase then becomes 



tp = 4> + <j) D = tp + Sip + <j) D 



(25) 



The normal to the detector plane z is at constant an- 
gle 6 Z = 7r/3 from the normal to the ecliptic z, and 
constantly points in the direction of the z-axis from the 
barycenter of LISA. Furthermore, each satellite rotates 
around the z-axis once a year (see Fig. [T|). Let us ex- 
press then the detector frame in the Solar System frame, 



FIG. 1: The orbit of LISA around the Sun as currently planned. Image taken from the LISA pre-phase-A report [2!|. 



assuming that y ■ y — 1 at t = 0: 
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LISA will act during the incoming of a gravitational 
wave as a pair of two-arm detectors, but with a response 
scaled by a a/3/2 factor due to the 60° opening angle of 
the constellation, following the pattern (k = 1,2): 
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where #at and </>jv are the spherical angles of the position 
of the binary in the detector frame, and ipN is defined 
through 
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We expressed here two combinations of the response of 
the three arms of LISA whose detector noises are uncor- 
rected H. 



Using this, we find the response function of the detec- 
tors (k = 1,2): 
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We can change this into the phase-amplitude represen- 
tation: 
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where is the polarization phase, and A?° n is the po- 
larization amplitude: 
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The final form of the gravitational wave signal is thus 
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To estimate the measurement error in the different 
parameters of the binary, we need to know the Fourier 
transform of the signal hk(f) [40( ] . 
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Note that we neglected in the last line the Fourier 
transform of the so-called memory effect A^?\ This is 
based on the fact that the Fourier transform of the func- 
tion xA+F£ accumulates around frequencies which are 
separated from the orbital frequency range, at least dur- 
ing most of the inspiral. It will thus not contribute to 
the relevant frequencies. 

To compute the integrals, we rely on the stationary 
phase approximation. Neglecting the integrals with the 
e i(2nft+ip kin ) f ac t orj as they will only contribute to nega- 
tive frequencies, the stationary points for the other inte- 
grals are given by 

2wf = 1p' k „(tfc,„) = 1Kj(t k , n ) + ra0£>(*k,n) + (<£fc° n )'(*fe,n)- 

(43) 

For the same reasons as before, we can safely neglect 
the derivatives of the Doppler phase and of the polariza- 
tion phase. We thus get the following expression for the 
stationary point: 



tk,n = tn = t(f/n), 



(44) 



where the function t(f) is defined at 2PN order by 
Eq. ©. 

Thus, we get the following expression for the Fourier 
transform of the gravitational wave signal: 
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where we used here x = x(u> = 2nf) a different orbital 
frequency parameter for each harmonic. 

Note that Lang and Hughes [l7| took the zeroth order 
form for S, S(f) = 1, consistently with neglecting all 
amplitude modulations. 

Finally, a binary will be observed with LISA during a 
finite amount of time. Therefore, if we denote by U and ir- 
respectively the initial and final time of observation, the 
orbital frequencies available for the Fourier transform will 
lie between f or b{ti) and forb(tf)- Thus, the final Fourier 
transform will be of the form: 

h{f) = E KnU)0{f - n/ort(ti))0(n/ort(t/) - /), 

where is the Heaviside step function. 

We compared in this work three different waveforms, 
which we called full waveform (FWF), simplified wave- 
form (SWF), and restricted waveform (RWF). The latter 
is the one used in (l7| . 

The FWF contains all post-Newtonian corrections of 
the frequency and amplitude of the wave up to 2PN or- 
der. It is obtained using the amplitudes given in the 
appendix in Eqs. (|38)) and (|39|) . and inserting the results 
in Eq. (g5J). 

The SWF contains all post-Newtonian corrections of 
the frequency up to 2PN order, and the lowest order 
amplitude of each harmonic present at the 2PN level. 
With this approximation, we find particularly simple 
forms for the polarization amplitudes and phases (with 
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The SWF is obtained inserting the polarization ampli- 
tudes and phases above into Eq. (|43|) and, consistently 
with neglecting all amplitude corrections, taking the low- 
est order of the overall amplitude correction S(f) = 1. 

The RWF contains all post-Newtonian corrections of 
the frequency up to 2PN order, and the lowest order am- 
plitude of the second harmonic. It is identical to the 
SWF, with the further approximation A%° = 4> p ^ l n = 
0, n ^ 2. "' " 



C. Data analysis 

The signal will of course also include noise. A good 
description of the impact of noise can be found in [17]], 
and we will refer to that study for how to model it. A 
deeper study of the Fisher information formalism in the 
context of gravitational wave experiments can be found 
in [H. 

As defined in [l7[, the inner product in the space of 
signals is 

W»)-«.jT ™* (50, 

where Sh(f) is the noise spectral density. 

Now, if we have a signal h, described by a certain set 
of n parameters 9, the Fisher information matrix is an 
n-by-n symmetric matrix defined by 
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When we have several detectors, the Fisher information 
matrices are simply added: 



r, 
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Finally, the covariance matrix is the inverse of the infor- 
mation matrix: 



e = r _1 . 



(53) 



Its off-diagonal elements represent correlation coeffi- 
cients between the different parameters, and must satisfy 



< 1, 
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whereas its diagonal elements represent lower limits to 
statistical errors on their measurement: 
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We chose to use a different noise curve than , moti- 
vated by the fact that an agreement seems to have been 
found among the LISA parameter estimation commu- 
nity [27| to use a piecewise fit of the expected noise, used 
in the second round of the mock LISA data challenge [3l| . 
We use the noise model described in (2?j], which consists 
in a sum of an instrumental noise S n (f) and a Galac- 
tic confusion noise S con {(f) 41]. The instrumental and 
confusion noise curves read (/ is given in Hertz) 



1 

1+ 2 



!_ 

/* 



•5 1) 



Sconf (/) 



^g-44.62 J-2.3 
^Q-50.92 J-4.4 
2Q-62.8 

lo- 89 - 68 /- 





4S a 

(27T/) 4 
-3\ 



(56a) 



-20 



(icr 3 < / < icr 2 - 7 ), 
(io- 2 - 7 < / io~ 2 - 4 ), 
(icr 2 - 4 < / 10- 2 ), 
(io- 2 < /). 

(56b) 



where L — 5 • 10 9 m is the arm length of LISA, S p = 

4 • 1CP 22 m 2 Hz -1 is the white position noise level, S a — 
9 ■ 10~ 30 m 2 s~ 4 Hz -1 is the white acceleration noise 
level, and /* = c/(2irL) is the arm transfer frequency. 
This way, we have S h (/) = S n (/) + S coni (f). 

We found substantial differences in the noise curve 
that we used with respect the one used by Lang and 
Hughes [l7|. We plot the ratio of the two curves in Fig. [2] 

Our noise is in general larger for frequencies below 
1CP 2 Hz, which corresponds to the maximum frequency 
of the RWF emitted by binaries with total mass M s» 

5 • 1O 5 M (f m ax is proportional to 1/M). Therefore we 
expect our results to be more pessimistic than the ones 
in [I7I ]. However, we do not expect the comparison be- 
tween the different waveforms to depend very much on 
the particular noise curve used in such a study. 

We also compared the two instrumental noises and the 
two confusion noises separately. The two confusion noises 
are similar below 10 -3 Hz and differ above this value. 
They become negligible compared to the instrumental 
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10 



10 



0.001 0.01 
/[Hz] 



0.1 



FIG. 2: Ratio of the noise curve we used over the one used by 
Lang and Hughes, as a function of the frequency. We added on 
top the total mass of a binary which emits its second harmonic 
at the corresponding maximum frequency. At high frequen- 
c i es (/ ^ 5 ■ 10 -2 Hz), the curve used by Lang and Hughes 
is inaccurate. There is a peak at / = 2 ■ 10~ 3 Hz, where our 
noise is 4 times larger, and which corresponds to the maxi- 
mum frequency of binaries with a total mass M ~ 2-lO 6 M . 
Between 10" 3 and 2 ■ 10" 4 Hz, we find good agreement, and 
for lower frequencies our noise is much larger. 



noise above 5 • 10 -3 Hz. The peak visible in Fig. [5] at 
/ = 2 • 10~ 3 Hz is due to the differences in confusion 
noises. The discrepancies come from the fact that af- 
ter the publication of [l7|, a better estimation of what 
fraction of low-mass binaries could be resolved and not 
contribute to the confusion noise has been made [32j and 
used in the noise model of (27[. The instrumental noise 
used in [ItJ is a low-frequency approximation based on 
the online sensitivity curve generator provided by S. Lar- 
son. The noise we used differs from the one used by 
Larson for frequencies below 10 -4 Hz. This comes from 
the fact that Larson assumes white acceleration noise, 
whereas the authors of [l?) additionally considered it to 
increase as 1/f below 10~ 4 Hz. 



III. SIMULATIONS 

As a set of 12 intrinsic plus 3 extrinsic parameters for 
our simulations, we used: 

(i) log 10 mi/M© and log 10 m-ijMQ , for the masses of 
the two black holes. 

(ii) fii = cosdi and 4>h f° r the spherical angles of the 
orbital angular momentum L at 7 = i. 

(iii) /ii = cos Q\ and <pi for the spherical angles of the 
spin of the first black hole Si at 7 = A. 

Xi — q^i \ Si\ f° r the dimensionless strength of 
the spin of the first black hole, which has to satisfy 
0<xi< 1- 



(v) (12 = cos 62, <p2, and xi 1 same for the second black 
hole as for the first one. 

(vi) t c , the time of coalescence. 

(vii) (fi c , the phase at coalescence. As this phase is ran- 
dom and its determination is not of any astrophys- 
ical interest, we can safely neglect constants in the 
orbital phase, in particular Sip from Eq. (|18[) . 

(viii) fi n — cos 9 n and <j) n , the spherical angles of the 
position of the binary in the sky. 

(ix) d,L , the luminosity distance between the source and 
the Solar System. 

All angles are taken in the frame tied to the distant 
stars. We fix the zero point of time by the beginning of 
the LISA mission. 

We performed Monte Carlo simulations where, given 
a set of parameters, we evolved the binary backwards in 
frequency starting from lo(j = 1/6) using a fourth order 
adaptive Runge-Kutta algorithm to find L(u>), Si(u>), 
S2(u>), and 5(p(u>). We stopped the simulations either at 
t = 0, or when the frequency of the highest harmonic 
had gone below the LISA band, for 6lo < 3 • 1CT 5 Hz. We 
chose to start at 7 = 1/6 because it is the radius of the in- 
nermost stable circular orbit (ISCO) for a Schwartzschild 
black hole. Of course, when dealing with a spinning black 
hole, the radius of the ISCO can vary between 7 = 1 
and 7 = 1/9, depending on the spin parameter of the 
black hole and on the orientation of the orbit, but the 
series may not converge when 7 is close to 1, so that we 
chose to consider the post-Newtonian expansion as ac- 
curate for 7 1/6 , which seems to be a good enough 
prescription [33j, [34j • 

We then put these functions inside the Fourier trans- 
forms of the waves (p4"5l) . 

Five derivatives dgihk out of the 15 needed can be 
found analytically. Three simple ones are: 



dt c 



2mfh k (6*,f), 



dh k (V,f) _ h k (6ij) 

= -i^nhk^lO 3 , /), 



dd L d L 



dif c 



(57) 
(58) 
(59) 



where hk, n is the n-th harmonic component of h k - The 
other two are the derivatives with respect to ji n and 4> n . 
The only quantities in Eq. (|45[) that depend on these 
parameters are A^ n , 4> P k ° n , Sip, and 

The derivatives which we could not find analytically 
have been computed numerically using the relation 

dh k (03,f) h k (W + f) - h k (0J - /) 



d8 l 



(60) 
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where e is a small displacement of the parameter 9 l . We 
used a constant value of e = 10~ 7 for every parameter, 
except for 4>i for which e was divided by 2 — 2|/^|, [ii 
(i G {1,2}) for which e was divided by 5\i, and <pi for 
which e was divided by 10x^(1 — |m»D- The formula is 
accurate up to 0(e 2 ). 

We computed the functions L(ui), Si(u>), S^w), and 
Sip(ui) for each displacement of the parameters. 

We then evaluated numerically the integrals 
(dgihk\dgihk) to find the Fisher information ma- 
trix. Each harmonic hk. n {f) is truncated if necessary 
to remain inside the LISA band, which we take to be 
[3 • 1CT 5 , 1] Hz. 

We added the contributions from both detectors, and 
then inverted the matrix numerically to find the statisti- 
cal error estimates. 

We found that in some extreme situations, when L ■ n 
gets close to 1, the Runge-Kutta method fails to converge 
when computing Sip, because 



dSip 
dw 



L ■ h 



1 



L ■ h 



L x n 



dL 



1 



L x n 



L n -> 1. 



(61) 



In those situations, we chose to compute 8(p(u)) whenever 
L ■ n is to close to 1 with an approximate value, which is 



6ip(uj + Suj) r 
5ip(uj) + angle 



yL(uj + Suj) x , (^L(u 



(62) 



We ran different sets of simulations fixing the rcdshift 
and the masses, and selected the other parameters ran- 
domly, using a flat distribution. The bounds to put for 
the Monte Carlo selection of the different parameters are 
clear, except for t c . We chose the following bounds, con- 
sistently with [l7[ : the lower bound for t c is for the phys- 
ical separation parameter 7 to be equal to 1/6 at t = 
(combining Eqs. ((2|) and (jSj), and the higher bound is 
t c = 3 yrs (this is in fact the minimum science require- 
ment of the mission), which we take to be the duration 
of the LISA mission, so that the coalescence is visible 
during it. We computed for each set the mean measure- 
ment errors for the parameters and signal-to-noise ratio 
(SNR), comparing the output for the RWF, the SWF, 
and the FWF defined at the end of Sec. IITbTI 



IV. RESULTS 

We ran different sets of simulations, each of them at a 
redshift of z = 1, varying the masses between O(1O 5 M ) 
and O(10 8 Mq), and the mass ratio between 1:1 and 1:10. 
We did not vary the redshift, because, as described in 
Sec. HH it cannot be measured by a gravitational wave 
experiment. Furthermore, with a redshift to luminosity 
distance relation fixed, the only parameter varying with 
the redshift for constant redshifted masses is the luminos- 
ity distance. The statistical errors in this case scale for 



all parameters as (1 -I- z)c?l, as this parameter appears 
only as an overall factor in the waveforms. For exam- 
ple, the statistical error estimates on the parameters for 
a system with mi = 5 • 10 6 M Q , m,2 = 5 • 10 5 M Q , and 
z = 1 are exactly the same as those for a system with 
mi = 2 ■ 1O 6 M , m 2 = 2 • 10 5 M Q , and z = 4 (same red- 
shifted masses), multiplied by 2.5g?l(z = 4)/dr,(z = 1). 
For binaries with masses higher than 10 7 M Q , the results 
for the RWF cannot be trusted, as the second harmonic 
spends typically only a few orbits inside the LISA band, 
and even no signal at all can be observed for 10 8 M Q bi- 
naries. Each of our sets of simulations consisted in over 
a thousand binaries. We performed a posteriori statisti- 
cal checks showing that the medians should be correctly 
estimated up to a few percent. 

We present here in tables for all samples and all in- 
teresting parameters a best-case measurement error (5% 
quantile), a typical error (the median), and a worst-case 
error (95% quantile). The parameters we are interested 
in are the (redshifted) individual masses of the black 
holes, shown in Tables U and |TT] their spin parameters, 
shown in Tables IIIII and IIV1 the principal axes of the lo- 
calization ellipse in the sky, shown in Tables [V] and IVIl 
and the (redshifted) luminosity distance to the source, 
shown in Table EU 

We followed [Tt} to present as sensible quantities for 
the sky localization the principal axes 2a and 2b of the 
ellipse enclosing the region outside of which there is a 1/e 
probability of finding the binary. 

For binaries for which no signal can be extracted from 
the data, we fixed the errors to 00. For errors apparently 
meaningless, such as Ad^/dL > 1 or 2a > 2ir, we still 
provide the error as computed, because it can give an 
indication on up to what redshift quantities can be com- 
puted using the scaling property of the error with respect 
to (1 + z)d L . 

We found that the binaries can roughly be separated 
into three classes: low unequal-mass binaries, low equal- 
mass binaries (M < 1O 7 M ), and high-mass binaries 
(M > 1O 7 M ). We discuss below these three distinct 
cases, and plot the estimated error distributions for a rep- 
resentative sample of each one of the three classes on the 
parameters Ami/mi, Axi, 2a, and Ad^/d^. The distri- 
butions for Ami) mi are similar to those for Ami /mi, 
those for A%2 to those for A%i , and those for 2b to those 
for 2a. 

In general, for lower-mass binaries and independently 
on the mass ratio, we find that the errors expected for 
extrinsic parameters using the FWF are ~ 1.5 times the 
ones expected for the RWF. This factor is ~ 1.2 compar- 
ing the SWF to the RWF. This is changed when consid- 
ering higher-mass binaries, because the second harmonic, 
the only one present in the RWF, spends very few cycles 
inside the LISA band. 

To discuss the mass limit above which no information 
can be extracted from a system anymore, we present at 
the end the proportion of systems for which the individ- 
ual masses and the luminosity distance can be measured 
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mi [M ] 


7712 [Mq] 


Ami /mi 






5% quantile 


Median 


95% quantile 






RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


3 ■ 10 5 


10 5 


1.85 ■ 1(T 4 


1.69 • 10~ 4 


1.44 ■ 10 -4 


8.06 • 10" 4 


6.71 • 10~ 4 


4.76 ■ 10~ 4 


8.17 ■ 10" 3 


2.54 ■ 10" 3 


1.58 • 10~ 3 


10 6 


10 5 


2.44 ■ 10~ 4 


1.99 ■ 10~ 4 


1.63 ■ 10" 4 


7.13 • 10" 4 


5.36 • 10~ 4 


4.25 ■ 10~ 4 


4.63 • 10" 3 


2.82 ■ 10" 3 


2.13 • 10~ 3 


10 6 


3 ■ 10 5 


4.08 ■ 10~ 4 


3.59 • 10~ 4 


2.90 • 10~ 4 


1.36 ■ 10~ 3 


1.10 • 10 -3 


8.01 ■ 10~ 4 


1.11 ■ 10~ 2 


3.95 ■ 10" 3 


2.60 • 10~ 3 


3 ■ 10 5 


3- 10 5 


1.69 ■ 10~ 4 


1.62 ■ 10~ 4 


1.19 ■ 10" 4 


1.24 • 10 -3 


1.17 • 10~ 3 


2.91 ■ 10~ 4 


1.14 ■ 10" 2 


8.82 ■ 10" 3 


6.52 ■ 10~ 4 


10 6 


10 6 


3.59 ■ 10 -4 


3.53 ■ 10~ 4 


2.53 ■ 10" 4 


2.48 • 10" 3 


2.42 • 10" 3 


7.16 ■ 10~ 4 


2.98 ■ 10" 2 


2.29 ■ 10" 2 


1.57- 10" 3 


10 7 


10 6 


1.21 ■ 10~ 3 


7.93 • 10" 4 


4.46 • 10" 4 


4.34 • 10" 3 


2.49 • 10" 3 


1.37 ■ 10~ 3 


4.21 ■ 10" 2 


1.12 ■ 10" 2 


5.09 • 10 -3 


10 7 


3- 10 6 


3.39 ■ 10~ 3 


1.89 • 10~ 3 


8.23 • 10~ 4 


1.59 ■ 10~ 2 


4.26 • 10" 3 


1.79 ■ 10" 3 


0.140 


1.18 ■ 10" 2 


5.86 • 10~ 3 


10 7 


10 7 


2.20 ■ 10~ 2 


1.21 ■ 10" 2 


2.04 ■ 10 -3 


0.213 


8.97- 10~ 2 


5.79 ■ 10" 3 


1.47 


0.825 


1.57- 10~ 2 


3 ■ 10 7 


10 7 


0.377 


8.64 ■ 10~ 3 


4.35 ■ 10~ 3 


1.01 


1.99 • 10~ 2 


9.48 ■ 10~ 3 


3.23 


5.35 • 10" 2 


2.43 • 10~ 2 


3 ■ 10 7 


3- 10 7 


3.74 


0.525 


5.55 ■ 10" 2 


23.1 


2.26 


0.120 


115 


9.05 


0.386 


10 s 


10 7 


oo 


9.67 ■ 10" 2 


9.45 ■ 10~ 2 


oo 


0.276 


0.246 


oo 


1.27 


1.00 


10 8 


3- 10 7 


oo 


0.896 


0.963 


oo 


2.57 


2.86 


oo 


53.6 


59.7 



TABLE I: Median, 5% and 95% quantiles of the estimated measurement errors on mi for different sets of binaries located at 
redshift z = 1, with low unequal-mass binaries on top, low equal-mass binaries in the middle, and high-mass binaries at the 
bottom. 



mi[M ] 


m 2 [M Q ] 


Am2/m2 






5% quantile 


Median 


95% quantile 






RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


3 ■ 10 5 


10 5 


1.50 ■ 10~ 4 


1.37- 10~ 4 


1.18 • 10" 4 


6.54 • 10" 4 


5.44 • 10~ 4 


3.87 ■ 10~ 4 


6.64 ■ 10" 3 


2.06 ■ 10" 3 


1.28 • 10~ 3 


10 6 


10 5 


1.76 ■ 10~ 4 


1.41 ■ 10~ 4 


1.18 ■ 10~ 4 


5.05 • 10~ 4 


3.78 ■ 10~ 4 


3.03 • 10~ 4 


3.27 ■ 10" 3 


2.00 ■ 10" 3 


1.50 ■ 10~ 3 


10 6 


3 ■ 10 s 


3.26 ■ 10 -4 


2.89 • 10 -4 


2.33 • 10" 4 


1.08 • 10" 3 


8.84 • 10~ 4 


6.41 • 10 -4 


8.91 ■ 10" 3 


3.15 ■ 10" 3 


2.09 • 10~ 3 


3 ■ 10 5 


3- 10 5 


1.64 ■ 10~ 4 


1.59 ■ 10~ 4 


1.21 ■ 10" 4 


1.25 • 10~ 3 


1.17 ■ 10 -3 


2.90 ■ 10 -4 


1.15 ■ 10" 2 


8.75 ■ 10~ 3 


6.55 • 10~ 4 


10 6 


10 6 


3.53 ■ 10" 4 


3.48 • 10~ 4 


2.56 • 10" 4 


2.52 • 10" 3 


2.41 • 10" 3 


7.13 • 10~ 4 


2.97 ■ 10" 2 


2.30 ■ 10" 2 


1.57- 10 -3 


10 7 


10 6 


1.19 ■ 10~ 3 


7.20 • 10~ 4 


4.25 • 10~ 4 


3.48 ■ 10~ 3 


1.84 • 10~ 3 


1.04 ■ 10~ 3 


2.94 ■ 10' 2 


7.89 ■ 10~ 3 


3.59 • 10~ 3 


10 7 


3- 10 6 


3.20 ■ 10" 3 


1.66 ■ 10~ 3 


7.45 ■ 10" 4 


1.33 • 10" 2 


3.56 • 10" 3 


1.54 • 10" 3 


0.110 


9.44 • 10" 3 


4.74 ■ 10~ 3 


10 7 


10 7 


2.20 ■ 10~ 2 


1.20 ■ 10~ 2 


2.05 ■ 10~ 3 


0.208 


9.15 • 10" 2 


5.79 ■ 10" 3 


1.49 


0.820 


1.57- 10~ 2 


3 ■ 10 7 


10 7 


0.412 


1.06 • 10~ 2 


5.78- 10" 3 


1.57 


2.34 • 10~ 2 


1.30 ■ 10~ 2 


4.91 


5.54 ■ 10" 2 


3.27 • 10~ 2 


3 ■ 10 7 


3- 10 7 


3.57 


0.557 


5.54 ■ 10" 2 


23.0 


2.21 


0.120 


108 


8.89 


0.386 


10 8 


10 7 


oo 


0.264 


0.336 


oo 


0.867 


1.05 


oo 


3.23 


3.95 


10 8 


3- 10 7 


oo 


3.19 


3.56 


oo 


9.75 


10.3 


oo 


145 


160 



TABLE II: Median, 5% and 95% quantiles of the estimated measurement errors on mi for different sets of binaries located at 
redshift z = 1, with low unequal-mass binaries on top, low equal-mass binaries in the middle, and high-mass binaries at the 
bottom. 



with 50% and 25% accuracy, for all samples. We also plot 
for different mass ratios the maximum redshift at which 
information can be extracted from a binary system, as a 
function of m\. 



We present then for each waveform how far the mea- 
surement of supermassive black hole mergers could help 
determining the Hubble diagram. To do so, we compute 
up to what redshift half of the systems can be localized 
inside the field of view of Hubble and/or XMM-Newton 
(see e.g. (35[) which we take to be 30' wide, with an error 
on g?l smaller than 10%. 



A. Low unequal-mass binaries 

We put in this class all systems with total mass smaller 
than 10 7 M Q , and with a mass ratio of at least 1:3. We 
chose to present as a representative sample systems with 
mi = 10 6 M o and m 2 = 3-lO 5 M . We plot the estimated 
distribution of the errors on m\ in Fig.[3l on xi i n Fig.|4l 
on the sky positioning in Fig. [5l and on di, in Fig. [S] 

For these systems, the gain in accuracy obtained in the 
determination of all interesting parameters with respect 
to the RWF is typically a factor ~ 1.5 for the FWF and 
a factor ~ 1.2 for the SWF. However, when the mass 
ratio is close to 1:3, the distribution of the errors on the 
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mi [M ] 


7712 [Mq] 


Axi 






5% quantile 


Median 


95% quantile 






RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


3 ■ 10 5 


10 5 


4.95 ■ 1(T 4 


4.11 ■ 10~ 4 


2.94 ■ 10 -4 


1.47 • 10~ 3 


1.17 ■ 10~ 3 


8.12 ■ 10~ 4 


9.63 ■ 10" 3 


4.83 ■ 10" 3 


3.32 • 10~ 3 


10 6 


10 5 


3.81 ■ 10~ 4 


2.72 ■ 10" 4 


1.97 ■ 10" 4 


8.98 • 10" 4 


6.11 • 10~ 4 


4.38 • 10 -4 


2.72 • 10" 3 


1.78 • 10" 3 


1.39 ■ 10~ 3 


10 6 


3 ■ 10 5 


8.68 ■ 10~ 4 


7.30 ■ 10~ 4 


4.97 ■ 10~ 4 


2.16 • 10~ 3 


1.69 • 10 -3 


1.20 ■ 10 -3 


1.05 ■ 10~ 2 


5.57 ■ 10" 3 


3.93 ■ 10~ 3 


3 ■ 10 5 


3- 10 5 


9.96 ■ 1(T 4 


9.62 ■ 10~ 4 


7.46 ■ 10" 4 


6.81 • 10~ 3 


6.45 • 10~ 3 


3.93 • 10" 3 


8.34 ■ 10~ 2 


7.08 ■ 10" 2 


3.75 ■ 10~ 2 


10 6 


10 6 


1.53 ■ 10" 3 


1.50 ■ 10~ 3 


1.17- 10" 3 


1.31 • 10" 2 


1.27 • 10~ 2 


6.82 • 10" 3 


0.218 


0.189 


7.51 ■ 10" 2 


10 7 


10 6 


1.37 ■ 1(T 3 


8.44 • 10~ 4 


4.58- 10" 4 


3.68 • 10" 3 


1.90 • 10" 3 


1.07 ■ 10" 3 


1.64 ■ 10" 2 


6.33 ■ 10 -3 


3.45 • 10" 3 


10 7 


3- 10 6 


3.87- 10~ 3 


2.41 • 10~ 3 


1.21 • 10~ 3 


1.47 • 10~ 2 


5.93 • 10" 3 


3.03 ■ 10" 3 


0.118 


2.15 ■ 10~ 2 


1.24 ■ 10~ 2 


10 7 


10 7 


0.108 


5.00 ■ 10~ 2 


2.04 ■ 10 -2 


1.20 


0.488 


0.136 


9.77 


5.45 


1.17 


3 ■ 10 7 


10 7 


0.438 


1.93 ■ 10~ 2 


9.36- 10~ 3 


1.87 


6.19 • 10~ 2 


3.31 • 10~ 2 


7.24 


0.269 


0.147 


3 ■ 10 7 


3- 10 7 


16.8 


2.04 


1.31 


83.6 


11.9 


5.35 


499 


61.4 


26.3 


10 s 


10 7 


oo 


0.256 


0.316 


oo 


1.10 


1.35 


oo 


4.65 


5.62 


10 8 


3- 10 7 


oo 


3.47 


4.19 


oo 


15.1 


16.3 


oo 


172 


185 



TABLE III: Median, 5% and 95% quantiles of the estimated measurement errors on xi f° r different sets of binaries located at 
redshift z — 1, with low unequal-mass binaries on top, low equal-mass binaries in the middle, and high-mass binaries at the 
bottom. 



mi[M ] 


m 2 [M Q ] 


AX2 






5% quantile 


Median 


95% quantile 






RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


3 ■ 10 5 


10 5 


8.05 ■ 10~ 4 


7.19 • 10" 4 


5.40 • 10~ 4 


3.23 • 10~ 3 


2.70 ■ 10~ 3 


1.91 ■ 10~ 3 


2.30 ■ 10" 2 


1.34 ■ 10" 2 


9.81 • 10~ 3 


10 6 


10 5 


1.71 ■ 10~ 3 


1.09 ■ 10~ 3 


8.34 ■ 10~ 4 


5.24 • 10 -3 


3.61 • 10~ 3 


2.78 ■ 10~ 3 


3.20 ■ 10" 2 


2.00 ■ 10" 2 


1.45 ■ 10~ 2 


10 6 


3 ■ 10 s 


1.41 ■ 10~ 3 


1.22 • 10~ 3 


9.01 • 10" 4 


4.72 • 10" 3 


3.91 • 10" 3 


2.78 ■ 10" 3 


2.65 ■ 10" 2 


1.63 ■ 10" 2 


1.20 ■ 10~ 2 


3 ■ 10 5 


3- 10 5 


1.02 ■ 10~ 3 


9.58 ■ 10~ 4 


7.83 • 10" 4 


6.66 ■ 10~ 3 


6.34 • 10~ 3 


3.95 • 10~ 3 


7.84 ■ 10" 2 


6.33 • 10" 2 


3.68 ■ 10" 2 


10 6 


10 6 


1.67 ■ 10~ 3 


1.66 • 10~ 3 


1.22 • 10~ 3 


1.27 • 10" 2 


1.25 • 10~ 2 


6.76 ■ 10~ 3 


0.217 


0.186 


7.07- 10 -2 


10 7 


10 6 


6.33 ■ 10~ 3 


3.68 • 10~ 3 


1.98 • 10~ 3 


3.31 • 10~ 2 


1.63 • 10~ 2 


9.66 ■ 10~ 3 


0.187 


7.57 ■ 10" 2 


4.01 • 10 -2 


10 7 


3 • 10 6 


7.35 ■ 10~ 3 


4.64 ■ 10~ 3 


2.42 ■ 10" 3 


3.32 • 10" 2 


1.72 • 10~ 2 


9.46 ■ 10" 3 


0.193 


6.87 ■ 10" 2 


4.13 • 10" 2 


10 7 


10 7 


0.111 


4.66 ■ 10~ 2 


1.94 ■ 10" 2 


1.22 


0.511 


0.130 


9.94 


5.60 


1.11 


3 ■ 10 7 


10 7 


0.594 


4.02 ■ 10 -2 


2.07- 10" 2 


4.44 


0.205 


0.104 


26.8 


0.960 


0.496 


3 ■ 10 7 


3- 10 7 


16.1 


2.20 


1.34 


83.0 


11.4 


5.43 


515 


59.5 


26.5 


10 8 


10 7 


oo 


1.07 


1.38 


oo 


11.7 


13.7 


oo 


52.8 


62.6 


10 8 


3- 10 7 


oo 


8.31 


9.25 


oo 


48.0 


50.4 


oo 


657 


700 



TABLE IV: Median, 5% and 95% quantiles of the estimated measurement errors on ^2 for different sets of binaries located at 
redshift z = 1, with low unequal-mass binaries on top, low equal-mass binaries in the middle, and high-mass binaries at the 
bottom. 



individual masses for the RWF has a relatively long tail 
of bad errors, which is absent for the SWF and FWF. 



The fact that including such extra structure as con- 
tained in the FWF fails to provide much extra accuracy 
can allow including extra parameters in the template. 
It has been recently suggested that the eccentricity of 
SMBH binaries could be significant in the last stages of 
the inspiral [36| . Thus, inserting eccentricity parame- 
ters [37| could be important. Furthermore, GW obser- 
vations could help constraining alternative gravity theo- 
ries @, S, @] . 



B. Low equal-mass binaries 

We put in this class all systems of equal-mass black 
holes, with total mass smaller than 10 7 M Q . We chose 
to present as a representative sample systems with mi = 
ni2 = 3 ■ 10 5 A/ Q . We plot the estimated distribution of 
the errors on mi in Fig. [7l on \i m Fig- El on the sky 
positioning in Fig. [SJ and on in Fig. [TTJJ 

In these cases, the errors on extrinsic parameters are, 
as for unequal-mass systems, improved by a factor ~ 1.5 
for the FWF with respect to the RWF. The errors on 
the spins are improved for the worst cases by a factor 
2-4, and typically by a factor 1.5 - 2 for the FWF with 
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mi[M ] 


mi [M Q ] 


2a ['] 






5% quantile 


Median 


95% quantile 






RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


3- 10 5 


10 5 


5.09 


4.96 


3.18 


20.2 


18.5 


12.7 


92.2 


85.8 


67.2 


10 6 


10 5 


8.67 


8.13 


7.14 


34.1 


28.2 


20.1 


124 


100 


82.5 


10 6 


3 ■ 10 5 


8.95 


8.61 


6.28 


31.4 


28.0 


19.8 


124 


110 


85.0 


3- 10 5 


3 ■ 10 5 


6.00 


5.81 


3.73 


26.6 


25.4 


18.6 


113 


109 


97.7 


10 6 


10 6 


8.45 


8.40 


5.42 


38.5 


37.6 


26.3 


158 


154 


129 


10 7 


10 6 


19.2 


15.5 


8.34 


64.2 


48.4 


23.2 


316 


199 


113 


10 7 


3 ■ 10 6 


19.5 


17.4 


7.31 


84.3 


65.2 


31.0 


461 


283 


158 


10 7 


10 7 


32.7 


27.4 


11.6 


202 


155 


77.7 


1360 


818 


496 


3- 10 7 


10 7 


169 


68.5 


24.4 


1500 


319 


133 


16600 


1550 


738 


3- 10 7 


3 ■ 10 7 


7910 


537 


363 


188000 


2590 


2570 


2890000 


14700 


16400 


10 s 


10 7 


oo 


998 


1300 


oo 


3380 


4400 


oo 


18200 


23800 


10 s 


3 ■ 10 7 


oo 


5900 


6510 


oo 


26300 


31000 


oo 


237000 


279000 



TABLE V: Median, 5% and 95% quantiles of the estimated measurement errors on the major axis of the localization ellipse in 
the sky for different sets of binaries located at redshift z — 1, with low unequal-mass binaries on top, low equal-mass binaries 
in the middle, and high-mass binaries at the bottom. 



mi[M ] 


m 2 [M ] 


26 ['] 






5% quantile 


Median 


95% quantile 






RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


3 ■ 10 5 


10 5 


0.795 


0.778 


0.453 


3.76 


3.49 


2.10 


13.8 


12.1 


7.40 


10 6 


10 5 


2.17 


1.69 


0.985 


10.2 


7.61 


4.50 


23.6 


15.5 


10.2 


10 6 


3- 10 5 


1.83 


1.63 


0.984 


8.63 


7.70 


4.47 


24.4 


19.2 


12.5 


3 ■ 10 5 


3 ■ 10 5 


1.00 


1.00 


0.575 


5.52 


5.29 


3.17 


20.3 


18.6 


13.8 


10 6 


10 6 


1.62 


1.61 


0.948 


9.11 


9.04 


5.26 


31.9 


29.9 


19.4 


10 7 


10 6 


3.40 


2.81 


1.17 


15.7 


12.5 


5.20 


41.3 


27.2 


12.3 


10 7 


3- 10 6 


3.04 


2.65 


1.11 


13.8 


11.9 


4.87 


54.3 


37.2 


17.3 


10 7 


10 7 


5.77 


4.64 


1.93 


24.8 


21.6 


9.49 


125 


90.8 


48.3 


3 ■ 10 7 


10 7 


36.4 


11.1 


4.00 


164 


47.9 


17.5 


896 


157 


69.9 


3 ■ 10 7 


3- 10 7 


727 


108 


57.2 


5740 


306 


230 


85000 


1090 


1350 


10 8 


10 7 


oo 


201 


251 


oo 


724 


930 


oo 


2610 


3600 


10 s 


3- 10 7 


oo 


1140 


1340 


oo 


3890 


4630 


oo 


44400 


52500 



TABLE VI: Median, 5% and 95% quantiles of the estimated measurement errors on the minor axis of the localization ellipse in 
the sky for different sets of binaries located at redshift z — 1, with low unequal-mass binaries on top, low equal-mass binaries 
in the middle, and high-mass binaries at the bottom. 



respect to the two other waveforms. However, the error 
on the individual masses is improved typically by a factor 
3.5 - 4.5, and even by a factor 10 - 20 in the worst cases, 
comparing the FWF with the two others. Thus, much 
more information can be extracted from a measure of 
an equal-mass binary system using the former waveform 
than one of the latter. 



The SWF brings little improvement for intrinsic pa- 
rameters in these cases, because the odd harmonics are 
absent from it, so that it has only two corrections to the 
RWF instead of five for unequal-mass systems. 



C. High- mass binaries 

We put in this class all systems with total mass higher 
than 10 7 Mq. We chose to present as a representative 
sample systems with mi = 3 • 10 7 M Q and = 10 7 M Q . 
We plot the estimated distribution of the errors on mi 
in Fig. 1111 on Xi m Fig- H21 on the sky positioning in 
Fig. [131 and on d L in Fig. QH 

For this class of binaries, the second harmonic is hardly 
or not at all visible in the LISA band, so that the RWF 
fails to provide good accuracy unless the total mass is 
close to 10 7 M Q . However, the SWF and FWF still 
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mi [M Q ] 


7712 [Mq] 


Ad L /d L 






5% quantile 


Median 


95% quantile 






RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


3 ■ 10 5 


10 5 


8.67 ■ 1(T 4 


8.24 • 10~ 4 


5.76 ■ 10~ 4 


2.94 • 10~ 3 


2.43 • 10~ 3 


1.75 ■ 10~ 3 


1.28 ■ 10~ 2 


1.01 ■ 10" 2 


8.43 • 10~ 3 


10 6 


10 5 


2.08 ■ 10~ 3 


1.44 ■ 10~ 3 


1.03 ■ 10" 3 


4.90 • 10" 3 


3.51 • 10" 3 


2.44 • 10" 3 


1.49 ■ 10" 2 


8.74 • 10" 3 


7.04 ■ 10~ 3 


10 6 


3 ■ 10 5 


1.75 ■ 10~ 3 


1.43 • 10 -3 


1.11 • 10~ 3 


4.62 ■ 10~ 3 


3.74 • 10 -3 


2.70 ■ 10 -3 


1.97 ■ 10~ 2 


1.39 ■ 10~ 2 


1.12 • 10~ 2 


3 ■ 10 5 


3- 10 5 


1.13 ■ 10~ 3 


1.08 ■ 10~ 3 


7.20 ■ 10" 4 


4.41 • 10~ 3 


3.80 • 10~ 3 


2.88 • 10~ 3 


1.89 ■ 10" 2 


1.48 • 10" 2 


1.27 ■ 10~ 2 


10 6 


10 6 


1.85 ■ 10~ 3 


1.77 ■ 10~ 3 


1.12 ■ 10" 3 


6.88 • 10" 3 


6.20 ■ 10" 3 


4.36 • 10" 3 


2.65 ■ 10" 2 


2.02 ■ 10" 2 


1.71 ■ 10" 2 


10 7 


10 6 


3.14 ■ 10~ 3 


2.32 • 10 -3 


1.48 • 10~ 3 


8.53 • 10" 3 


6.25 • 10" 3 


3.45 • 10~ 3 


4.14 ■ 10" 2 


1.94 ■ 10" 2 


1.26 ■ 10~ 2 


10 7 


3- 10 6 


3.84 ■ 10~ 3 


3.11 • 10~ 3 


2.18 • 10~ 3 


1.09 ■ 10~ 2 


8.03 • 10" 3 


4.69 ■ 10" 3 


7.26 ■ 10~ 2 


3.43 ■ 10" 2 


1.99 ■ 10~ 2 


10 7 


10 7 


1.45 ■ 10~ 2 


8.99 ■ 10~ 3 


6.78 ■ 10" 3 


6.05 • 10" 2 


2.76 ■ 10~ 2 


2.03 • 10~ 2 


0.277 


0.113 


7.82 ■ 10" 2 


3 ■ 10 7 


10 7 


0.212 


1.94 ■ 10~ 2 


1.68 ■ 10~ 2 


0.801 


4.63 • 10~ 2 


3.52 • 10~ 2 


3.44 


0.186 


0.101 


3 ■ 10 7 


3- 10 7 


3.52 


0.188 


0.213 


31.7 


0.614 


0.523 


462 


2.33 


2.23 


10 s 


10 7 


oo 


0.257 


0.424 


oo 


0.698 


1.11 


oo 


2.53 


4.00 


10 8 


3- 10 7 


oo 


1.76 


3.17 


oo 


5.53 


9.16 


oo 


69.7 


111 



TABLE VII: Median, 5% and 95% quantiles of the estimated measurement errors on cLl for different sets of binaries located 
at redshift 2 = 1, with low unequal-mass binaries on top, low equal-mass binaries in the middle, and high-mass binaries at the 
bottom. 
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RWF 
SWF 
FWF 



0.01 



0.1 




FIG. 3: Estimated distribution of the measurement error on 
mi for a low unequal-mass binary system with mi = 1O 6 M0 
and m2 = 3 ■ 1O 5 M0. We expect to have errors as high as 
8 ■ 10 -3 with the RWF (95% percentile), whereas we do not 
expect errors higher than 1.5 • 10 -3 with the FWF. 

provide relatively high precision measurement on the 
masses, spins and luminosity distance of a system with 
mi — 3 ■ 1O 7 M at redshift one. 

For equal-mass systems in this class, the FWF pro- 
vides in all cases an improvement of a factor 10 - 30 
for the determination of the masses with respect to the 
SWF. For other parameters and/or other mass ratios, the 
improvement using the FWF is of a factor 1.5-2 with 
respect to the SWF. The fact that the SWF seems to 
give better results than the FWF for the highest-mass 
systems comes from the fact that the SNR for systems 
in this mass range is higher with the former than with 
the latter. However, we do not expect to extract more 



FIG. 4: Estimated distribution of the measurement error on 
Xi for a low unequal-mass binary system with mi = 1O 6 M0 
and ni2 = 3 • 10 Mq. We expect the error to be 1.5 - 2 times 
lower using the FWF than using the RWF. 



information from a more approximate waveform. 

Furthermore, some information can still be extracted 
from binaries that are completely invisible to the RWF 
using higher harmonics. 



D. Upper mass limit 

We present here for all samples, the proportion of sys- 
tems for which both individual masses can be measured 
at the level of 25% and 50% at a redshift of z = 1 in Ta- 
ble Ivnil as well as the luminosity distance in Table Hxl 
When at least 25% accuracy is obtainable for all systems 
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FIG. 5: Estimated distribution of the major axis of the po- 
sitioning error ellipse for a low unequal-mass binary system 
with mi = 10 6 Mq and m.2 = 3 • 10 5 M©. We expect the error 
to be ~ 1.5 times lower using the FWF than using the RWF. 



FIG. 7: Estimated distribution of the measurement error on 
mi for a low equal-mass binary system with mi = m,2 = 
3- 1O 5 M . The FWF clearly gives better results than the two 
other waveforms. 
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FIG. 6: Estimated distribution of the measurement error on 
dh for a low unequal-mass binary systems with mi = 10 6 Mq 
and ni2 = 3-10 M©. We expect the error to be ~ 1.5 times 
lower using the FWF than using the RWF. 



FIG. 8: Estimated distribution of the measurement error on 
Xi for a low equal-mass binary system with mi = m,2 = 
3 • 10 Mq. The improvement on the median value is a factor 
1.5 - 2 and on the 95% quantile a factor 2 - 4 for the FWF 
with respect to the other two waveforms. 



of a sample, we do not show it on the table. 



mi [M©] 


m 2 [M ] 


25% 


50% 






RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


10 7 


3- 10 6 


98% 


100% 


100% 


100% 


100% 


100% 


10 7 


10 7 


50% 


75% 


100% 


71% 


89% 


100% 


3- 10 7 


10 7 


1% 


100% 


100% 


5% 


100% 


100% 


3- 10 7 


3- 10 7 


0% 


0% 


84% 


0% 


3% 


98% 


10 8 


10 7 


0% 


4% 


2% 


0% 


21% 


15% 


10 8 


3- 10 7 


0% 


0% 


0% 


0% 


0% 


0% 



TABLE VIII: Proportion of the systems in all samples where 
both individual masses can be determined with an accuracy 
better than 25%, resp. 50%. 



mi[M ] 


m 2 [M Q ] 


25% 


50% 






RWF 


SWF 


FWF 


RWF 


SWF 


FWF 


10 7 


10 7 


94% 


99% 


100% 


99% 


100% 


100% 


3- 10 7 


10 7 


8% 


97% 


99% 


28% 


100% 


100% 


3- 10 7 


3 ■ 10 7 


0% 


12% 


9% 


0% 


40% 


47% 


10 s 


10 7 


0% 


5% 


1% 


0% 


29% 


10% 


10 8 


3 ■ 10 7 


0% 


0% 


0% 


0% 


2% 


0% 



TABLE IX: Proportion of the systems in all samples where 
the luminosity distance can be determined with an accuracy 
better than 25%, resp. 50%. 



We see in the tables that the RWF reaches its limits 
for 10 7 M Q binaries, whereas the SWF and FWF can still 
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FIG. 9: Estimated distribution of the major axis of the posi- 
tioning error ellipse for a low equal-mass binary system with 
mi = m2 = 3 ■ lO 5 A/0. We expect the error to be ~ 1.5 times 
lower using the FWF than using the RWF. 
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FIG. 10: Estimated distribution of the measurement error 
on (1l for a low equal-mass binary system with mi = 777,2 = 
3 • 10 Mq. We expect the error to be ~ 1.5 times lower using 
the FWF than using the RWF. 



provide significant information for 3T0 7 M Q binaries, and 
even for some 1O 8 M0 binaries with high enough mass 
ratio. 

Furthermore, we computed from our simulations the 
maximum redshift at which a binary system is observ- 
able, as a function of m l5 for different values of the 
mass ratio. We chose to call a system with parameters 
(mi, m2,;z) observable if at least half of the systems of 
these masses at this redshift have both individual masses 
measurable with at least 25% precision. We present in 
Fig. [15] the maximum redshift at which equal-mass sys- 
tems can be observed, in Fig. [T75] the same for systems 
with a mass ratio between 1:3 and 3:10, and in Fig. [T7] 
the same for systems with a mass ratio of 1:10. Some 
points are absent of the plots, because no signal at all 
could be extracted from the RWF when the higher-mass 



FIG. 11: Estimated distribution of the measurement error on 
mi for a high-mass binary system with mi = 3 ■ 1O 7 M0 and 
7772 = 1O 7 M0. Very few systems are measurable with the 
RWF with a precision less than 50%, whereas the other two 
waveforms provide in the worst cases a few percent precision, 
the FWF typically a factor of 1.5 - 2 better than the SWF. 




0.001 



FIG. 12: Estimated distribution of the measurement error 
on xi f° r a high-mass binary system with mi = 3 • 1O 7 M0 
and 777.2 = lO 7 Af0. We can see that no information on the 
spins can be extracted with the RWF, whereas some can be 
extracted with the two others in all cases, a factor of two 
better for the FWF than for the SWF. 



black hole had a redshifted mass mi » 1O 8 M0. 

The figures show that a much higher redshift can be 
reached with the FWF than with the other waveforms, 
and that the difference is bigger for mass ratios closer to 
1:1. The FWF is giving the possibility to observe any 
binary system with total mass of M ^ 10' Mq up to a 
redshift of z > 10, whereas the other waveforms fail, es- 
pecially for equal-mass systems. The same combinations 
of redshifted masses can be observed with the FWF at 
redshifts 1.5-5 times higher than with the RWF, which 
could greatly help constraining black hole and galaxy for- 
mation models [l|. 
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FIG. 13: Estimated distribution of the major axis of the 
positioning error ellipse for a high-mass binary system with 
mi = 3 ■ 1O 7 M0 and 7712 = 10 7 Mq. We expect to have a 
positioning error in the best cases (5% quantile) of 2.8° for 
the RWF, of 1° for the SWF, and of 25' for the FWF. 




FIG. 14: Estimated distribution of the measurement error 
on d,L for a high-mass binary system with mi = 3 • 1O 7 M0 
and 7712 = 10 7 A/q. We do not expect a measurement more 
accurate than 20% to be possible with the RWF, whereas the 
accuracy should be always less than 20% with the SWF, and 
less than 10% with the FWF. 



E. Extrinsic parameters 

We plot here as a function of the mass of the most 
massive black hole the maximum redshift where the ma- 
jor axis of the positioning error ellipse is less than 30' for 
half of the binaries. Equal-mass binaries are shown in 
Fig. 1181 binaries with mass ratio between 1:3 and 3:10 in 
Fig. [THl and binaries with a mass ratio of 1:10 in Fig. [201 

We found that in all cases, the localization in the sky- 
is far more difficult than the determination of the lu- 
minosity distance. Irrespective of the waveforms, masses 
and mass ratios, the luminosity distance can be measured 
with a precision of 0.3% - 0.5% when the major axis of 



FIG. 15: Maximum redshift at which a system is observable 
as a function of the mass of the black holes, for equal-mass 
systems. The FWF allows to observe 10 4 - 1O 5 M0 binaries 
up to z = 30 - 50, the other two up to 2 — 15 - 25. A binary 
of ~ 2 • 10 Mq black holes should be observable up to z ~ 2 
with the FWF, awl with the SWF, and 2 w 0.3 with the 
RWF. 



z 




0.1 1 — 1 — 1 — 1 —J — 1 "J — 1 ^ — 1 —J 
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rrq [M ] 

FIG. 16: Maximum redshift at which a system is observable 
as a function of the mass of the most massive black hole, for 
systems with a mass ratio between 1:3 and 3:10. A binary 
with mi « 2 • 1O 4 M0 should be observable up to 2 ~ 36 with 
the FWF, uptozsi 30 with the SWF, uptozRi 28 with the 
RWF. A binary with mi w 1O 7 M0 should be observable up 
to z « 7 with the FWF, up to z w 6 with the SWF, and up 
to 2 w 2 with the RWF. 

the positioning error ellipse is 30'. 

The FWF could help locating binaries accurately 
enough for the observation of their merger to become 
possible up to redshifts 0.2 - 0.4 greater than the two 
other waveforms. The SWF could furthermore, in the 
case of unequal- mass binaries, go up to redshifts ~ 0.1 
greater than the RWF. 

Our simulations show that supermassive black hole bi- 
naries could be very accurate standard candles, and could 
successfully extend the measurements of the Hubble di- 
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FIG. 17: Maximum redshift at which a system is observable 
as a function of the mass of the most massive black hole, 
for systems with a mass ratio of 1:10. A binary with mi ~ 
6 • 1O 4 M should be observable up to z « 38 with the FWF, 
up to z « 34 with the SWF, up to z « 29 with the RWF. 
A binary with mi ~ 10 8 Mq should still be observable up to 
z « 0.7 with the FWF and SWF, and not visible at all with 
the RWF. 




FIG. 18: Maximum redshift at which the binary can be lo- 
cated with a 30' precision as function of the mass of the most 
massive black hole, for equal-mass systems. The FWF allows 
to locate a binary accurately up to a redshift of ~ 0.2 greater 
than the two other waveforms. 



agram up to rcdshifts of z = 1.6, with a precision on the 
luminosity distance of a few per mille. This would be 
a great breakthrough in the distance ladder, as the cur- 
rent most effective standard candles at large distances, 
type la supernovae, provide much less precision on large 
luminosity distances. 



FIG. 19: Maximum redshift at which the binary can be lo- 
cated with a 30' precision as a function of the mass of the most 
massive black hole, for systems with a mass ratio between 1:3 
and 3:10. The FWF allows to locate a binary accurately up 
to a redshift of 0.2 - 0.3 greater than the SWF, and the SWF 
up to a redshift of ~ 0.1 greater than the RWF. 
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FIG. 20: Maximum redshift at which the binary can be lo- 
cated with a 30' precision as a function of the mass of the 
most massive black hole, for systems with a mass ratio of 
1:10. The FWF allows to locate a binary accurately up to a 
redshift of 0.2 - 0.4 greater than the SWF, and the SWF up 
to a redshift of ~ 0.1 greater than the RWF. 



V. CONCLUSION 

The fact that the detection of gravitational waves 
with interferometric detectors relies on template-based 
searches should suggest to use the most accurate wave- 
form available for detecting systems emitting such waves. 
The gravitational waveform of two spinning bodies orbit- 
ing each other is however complicated, and each new fur- 
ther step implies more and more complicated corrections 
to the waveform. It is a good thing to know whether it 
is worth using a more accurate waveform, and in what 
cases. 
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Comparing our results for the RWF to those of Lang 
and Hughes [lj} , we found that our estimates are system- 
atically a factor of ~ 1.2 more pessimistic for lower- mass 
binaries, and a factor of ~ 2 more pessimistic for higher- 
mass binaries. As stated in section III C[ we expected 
such a discrepancy because of the differences in the noise 
curves we used. As higher-mass binaries spend more time 
in the lower frequency range where the noise differs the 
most, we expected the differences for these binaries to be 
larger than for less massive binaries. 

We found that the addition of higher harmonics to the 
waveform at the 2PN level can help increasing the mass 
limit above which no information can be extracted from 
the signal. The amplitude corrections can also bring im- 
portant improvements to the determination of the indi- 
vidual masses, also for lower-mass binaries. The FWF al- 
lows for detecting binaries up to redshifts z > 40, whereas 
the other waveforms can allow detection up to z ~ 30. 
This could be very interesting for constraining galaxy and 
black hole formation models. 

The range of LISA could be also extended for the de- 
termination of the Hubble diagram. The RWF would al- 
low to measure the redshift-luminosity distance relation 
at a few per mille precision up to z ~ 1.2, whereas the 
FWF would allow measures with the same precision up 
toz« 1.6. It would be interesting to quantify how well 
LISA would be able to determine the Hubble diagram. 

The use of the full waveform as a template for the 
gravitational radiation of comparable-mass binaries can 
be important to extract the maximum information pos- 
sible, especially for high-mass and/or close to equal-mass 
binaries. However, in the case of unequal low-mass bi- 
naries at low redshifts, the restricted waveform used in 
earlier studies can be sufficient. The fact that using the 
full waveform in these searches can fail to provide much 
more accuracy for some systems suggests that including 
more parameters, such as eccentricity [37| or alternative 
gravity parameters [3, S II] > could keep the accuracy for 
the other parameters reasonably high, allowing to extract 
more information from the detection of a wave. Gravi- 
tational waves can be a powerful tool for constraining 
alternative theories of gravitation, in the sense that each 
event will provide an independent measurement of their 
parameters. 

Even though the spin-coupling effects in the wave am- 
plitude are not yet known at the 2.5PN level, it could 
be interesting to compare the measurement accuracy we 



get from a 2.5PN accurate waveform as compared to the 
2PN accurate one we used in this study. It has been 
shown [H, [3§| in the case of spinless bodies that theoret- 
ical errors due to the inaccuracy of the waveform can be 
important for high SNR systems. It would also be inter- 
esting to perform the same study for spinning systems. 
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APPENDIX: POLARIZATIONS 

We give here the plus and cross polarizations we used 
in our studies, in terms of the orbital phase ip. The plus 
and cross polarizations of the simplified waveform (SWF) 

is obtained by keeping only the lowest order in and 

(n) 

B x , and those of the RWF by keeping only the lowest 

order of and . To obtain the SWF and RWF, the 
function S(f) in Eq. (05} should also be set to S(f) = 1. 
The plus and cross polarization waveforms are: 
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With the use of the spin-orbit coupling parameter r 
defined as: 
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